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Abstract 

We use the theory of the fluctuating electromagnetic field to cal- 
culate the frictional drag between nearby two-and three dimensional 
electron systems. The frictional drag results from coupling via a fluc- 
tuating electromagnetic field, and can be considered as the dissipative 
part of the van der Waals interaction. In comparison with other simi- 
lar calculations for semiconductor two-dimensional system we include 
retardation effects. We consider the dependence of the frictional drag 
force on the temperature T, electron density and separation d. We 
find, that retardation effects become dominating factor for high elec- 
tron densities, corresponding thing metallic film , and suggest a new 
experiment to test the theory. The relation between friction and heat 
transfer is also briefly commented on. 

1 Introduction 

A great deal of attention has been recently devoted to double layer systems in 
which two parallel quasi-two-dimensional (2D) subsystems (electron or hole 
gases) are separated by a potential barrier thick enough to prevent particles 



from tunneling across it but allowing for the interaction between particles on 
both its sides. Some time ago, Pogrebinskii and later Price 

predicted that the Coulomb interaction between two 2D electron sys- 
tems will induce a frictional drag force between the layers: a current in one 
film will induce a current in the adjacent film. The first frictional drag ex- 
periment was performed by Gramila et al. for two electron layers [21 El and 
by Sivan, Solomon and Shtrikman for an electron- hole system j^J. In these 
experiments a current is drawn in the first layer, while the second layer is 
an open circuit. Thus no dc-current can flow in the second layer, but an 
induced electric field occur that opposes the " drag force " from the first 
layer. These experiments spurred a large body of theoretical work both on 
electron-hole systems |5] and on electron-electron systems 0-113 • Most of 
this work focused on interlayer Coulomb interaction, the most obvious cou- 
pling mechanism and the one considered in the original theoretical papers 
PP, though the contributions due to an exchange of phonons between the 
layers have also been considered El M EE] ■ The origin of Coulomb drag is 
due to quantum and thermal fluctuations of the charge and current densities 
and can be considered as the dissipative part of van der Waals interaction. 
The static aspects of van der Waals interaction is well understood, and from 
the theory of Lifshitz J7j it is known that one must distinguish between two 
distance regime : (a) The nonretarded limit, when the separation between 
bodies d is small compared to wavelengths A ~ c/u , where u> is a char- 
acteristic frequency of the charge fluctuation, and c the light velocity, the 
interaction is determined by the fluctuations in an instantaneous Coulomb 
field . For metal u ~ u p , where u p is the plasma frequency .(b) Retardation 
effects become important, when d > A . As we have shown in Ref. ^Hl 
when calculating the dissipative part of the van der Waals interaction for 
two semi-infinite bodies in relative motion, retardation effects become im- 
portant for d > c/ujq, where ojq ~ uj p (uj p t) and r is the relaxation time. For 
Up ~ 10 16 s _1 and r ~ 10~ 14 s retardation effects become important for very 
short distances d > 1 A. However for 2D systems there is no investigation 
of the role of retardation effects in the frictional drag experiments. For large 
distances the retarded contribution to the frictional drag becomes important, 
and it is interesting to compare this contribution to the non-retarded contri- 
bution. To evaluate the retarded contribution from photon exchange we use 
the general theory of the fluctuating electromagnetic field developed by Ry- 
tov [20] and applied by Lifshitz 17 for studying the conservative part, and by 
us [TH] for studying dissipative part of the van der Waals interaction. In this 
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approach the interaction between the bodies is mediated by the fluctuating 
electromagnetic field which is always present in the vicinity of any collection 
of atoms. Beyond the boundaries of a solid this field consist partly of trav- 
eling waves and partly of evanescent waves which are damped exponentially 
with the distance away from the surface of the body. The method we use 
for calculating the frictional drag between two nearby 2D systems is quite 
general, and is applicable to any body at arbitrary temperature. It takes 
into account retardation effects, which become important for large enough 
separation between the bodies. 

We shall calculate frictional stress a — acting on the electrons in 
layer 1 due to the current density J 2 = n 2 ev in the layer 2, where n 2 is 
the carrier concentration (per unit area). If no current is allowed to flow in 
layer 1 (open circuit) an electric field E% develops whose influence cancels 
the frictional stress a between the layers. The frictional stress o = 'jv must 
equal the induced stress ri\eE\ so that 

7 = n 1 eE 1 /v = n 1 n 2 e 2 E 1 / J 2 = riin 2 e 2 p l2 , 

where the transresistivity p 12 = E\j ' J 2 is defined as the ratio of the induced 
electric field in the first layer to the driving current density in the second 
layer. The transresistivity is often interpreted in terms of a drag rate which, 
in analogy with a Drude model, is defined by r^ 1 = p 12 n 2 e 2 /m* = cr/n 1 m*t>. 

We find that for modulation-doped semiconductor quantum wells retarda- 
tion effects are not important under typical experimental conditions, support- 
ing earlier calculations where retardation effects always have been neglected 
0, jol-fTS] . However, although previous calculation for friction drag be- 
tween two-dimensional semiconductor systems are equivalent to ours, other 
approaches were very different. The present derivation offers an alterna- 
tive insight and is more general. A striking new result we find, that for 
systems with high 2D-electron density, e.g., thin metallic films, retardation 
effects becomes crucial and in fact dominates the frictional shear stress a. 
To test the theoretical predictions presented below we therefore suggest per- 
forming experiments on thin metallic layers grown on insulating substrates 
and separated by thin insulating layers. For example, for two thin (of order 
monolayer) silver films separated by d ~ 100 A, we estimate that the induced 
voltage Ui in metal film 1, due to a current J 2 in layer 2, will be of order 
U\ ~ 10 -8 U 2 , where U 2 is the driving voltage applied to metal film 2. Thus if 
U 2 ~ IV the induced voltage will be of order 10 pV which should be possible 
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to detect experimentally. We note that the study of this problem is also of 
direct interest in the context of sliding friction, since the electronic friction 
probed when two metallic bodies slide relative to each other, should be the 
same as the electronic friction probed by the transresistivity measurement, 
see Fig. 1. The electronic sliding friction (usually called vacuum friction) 
has recently been invoked to explain experimental results for the damping of 
a small metal particle vibrating in the vicinity of a flat metal surface , 
but this explanation is controversial jTHj , and it is clear that independent 
studies of the electronic friction would be of great interest. 

2 Calculation of the fluctuating electromag- 
netic field 

Let us firstly calculate the fluctuating electromagnetic field from one 2D 
system, surrounded by a dispersionless dielectric medium. We introduce a 
coordinate system with the xy— plane in the 2D layer. Following Lifshitz 
[17J, to calculate the fluctuating field we shall use the general theory due 
to Rytov, which is described in his book [20] • This method is based on 
the introduction of a "random " field in the Maxwell equations (just as, for 
example, one introduce a "random" force in the theory of Brownian motion of 
a particle ).For a monochromatic field (time factor exp(— iut)) in a dielectric, 
nonmagnetic medium, these equations are: 

V x E = i-B 

c 

VxH = -i-D+— (j+i f )5(z) (1) 

c c J 

where, following to Rytov, we divided the total current density j tot into two 
parts, jtot — j+j/i the fluctuating current density j/ associated with thermal 
and quantum fluctuations, and the current density j induced by the electric 
field E. D, H and B are the electric displacement field, the magnetic and 
the magnetic induction fields, respectively. For nonmagnetic medium B = H 
and D = eE, where e is the dielectric constant of the surrounded media. 
Eliminating B and H from one get 

V 2 E + (^)%E - ^V(p + Pf )5(z) + ^(j + j f )6(z) = (2) 
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where the total charge density p+Pf is the sum of the induced and fluctuating 
electron densities. We represent the current and electron densities in the form 
of a Fourier series 

j( r ) = 4rEj(qy q - r (3) 



A q 

p(r) = ^Ep(qK q - r (4) 

where q and r are 2D vectors in the xy-pl&ne and A is the surface area. From 
the equation of continuity 

one get up = q • j. Then (J2J) takes the form 

d 2 E(q) 2 Aniq . Attz . , 

+p 2 E q -q-(j +J/)<*(*) q- J +J/ )S'(z) 

d z sou 1 euj 1 



+ ^(i+j/W*) = o (5) 



C^ 2 



where 



and z is a unit vector along the z -axis. We shall consider separately the 
two cases where the electric field E is in the plane determined by the vectors 
q = q/g and z (p— polarized waves) and perpendicular to this plane along 
the vector n = z x q (s— polarized waves). 

Let us firstly suppose that E is parallel the vector n. In this case (J5|) 
gives 

^ + P 2 ^n + — ^ (<r t (q, c)£ n + j fn ) 5(z) = (7) 

where we have used Ohm's law 

in = cr t (q,uj)E a (8) 

where at is the transverse conductivity of the layer. The solution of (J7J) can 
be written in the form 

E n = u n e ip W (9) 
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^From (J7J) one can get the following boundary conditions at z = 



L 



E n {z = +0) = E n (z = -0) 

d^n d^n , 4njUJ f X 

— |,=+o - — | z =-o = — (a t (q,^)£ n + j /n J (10) 

Substituting into boundary conditions (fTUj) we get the following relation 
between the fluctuating current density jf n and electric field E n 

Jfn = -(cr t + ^-)E n (z = 0) (11) 

On the other hand the electric field can be calculated using linear response 
theory [221 • The Hamiltonian of the system has the form H + Hi nt , where H 
is the Hamiltonian of the body and radiation field, while 

H in t = -- / i n (r)j /n (r)e— l d 2 v = - ~ £ i n (-q)j/ n (q^K^ (12) 

C J c q 

Accordingly to the linear response theory the average value of the vector 
potential is determined by 



A n (q,w)) = at(q,uj)jf n (q^)/c (13) 

Taking into account that E = (iu/c)A, from comparison (jllj) and ()13)) we 
obtain 



Using fluctuation-dissipation theorem [23] we get 

+ n(u;)J Recr t 5 qq /5(co' — u'), (15) 



V V2 

where the Bose-Einstein factor 

n ^ = 6 ^a!t_i ' 

T is the temperature and Rea is the real part of conductivity. For q < 
{uj/c) yfe (fTHjl includes a term which does not depend on the conductivity, 
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and which result from the second term in the expression (|14j) . This term is 
associated with the black-body radiation which exist in the dielectric even 
without the 2D layer. Thus this term is irrelevant to the problem under 
consideration, and must be omitted. For p— polarized waves we get a similar 
expression as (fTo]) . except we must replace a t — > ui, where o\ is the longi- 
tudinal conductivity of the layer. Since the s— and p— polarized waves are 
not coupled, the average value of the product ^j/ q (q, c<j)j^ n ,(q / , u>')^ = 0. 

Let us now consider two parallel 2D electron layers separated by a distance 
d. We introduce two reference systems K and K', with coordinate axes xyz 
and x'y'z' . The xy- and x'y'- planes coincide with layer 1, with the x- and 
x' - axes pointing in the same direction, and the z- and z'- axes pointing 
toward layer 2. In the K system both layers are at rest. Assume now 
that in the layer 2 the conduction electrons move with the drift velocity 
v, corresponding to the current density ji = n,2ev, while no current flow in 
layer 1. The K' reference system moves with velocity v along to the x— axis 
relative to frame K. In the K' frame there is no current density in layer 2, 
while the surrounding dielectric moves with velocity —vx. In the K frame 
for z < d the Maxwell equations have the form with j and replaced 
by ji and respectively. After decomposition of the components of the 
electromagnetic field into a Fourier series the general solution of the Maxwell 
equation for z < d can be written in the form 

f ve ipz + we~ ipz , < z < d 
= u ie -^, z < (16) 



B 



( [q x v] + p [z x v])e ip2 + ([q x w] — p [z x w])e ipz , < z < d 
([qxu[]-j)[zx U!])e~ ip2 , z < 

(17) 

where v, w and ui satisfy the transversality conditions 

v ■ q+pv z = 0, w • q— pw z = 0, u r q+pu l2 , = (18) 

We now decompose the electromagnetic field into s— and p— polarized waves. 
The boundary conditions at z = for s— polarized waves is determined by 
(jiup. For p— polarized waves, from (JJJ) one obtain the boundary conditions 

E^z = +0) = E^z = -0) 
d£ q , d£ q , 4vrip 2 



, z=+0 . 



a 1 i(q,uj)E cl + j n(l ) (19) 
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^From (|1U|) and ()19|) we can obtain the following equations: 



w n + i?i s w; n = — — J ——r (21) 

pc 2 {e ls + 1) 

where v q = q • v and so on, and 

_ ei,(p) - 1 _ 4iru<7t _ inpcri 

eia(p) + 1 PC 2 we 

The Maxwell equations in the K' — system for z > have the same form 
as with j — > j 2 and jj — > jy 2 - However, to first order in v/c the relations 
between D, E, and B, H are |24] 



D=eE-(e-l)^xxB (22) 

H = B-(e-l)-xxE (23) 
c 

After eliminating D, B and H from Maxwell equations, and writing E, j and 
jf in Fourier series we get 

d 2 E'(q') /2 47rij/ 2 „, 4vrz , 

+P E q — q • j 2 +J /2 e q .<5 (z-d) q • J 2 +J /2 )$ {z-d) 

d -2 ecu ecu J 



+ ^^2+$^ = 

<r ec 

(24) 

where 



u '\ 2 '2 / I 2/5g x c 
- £ — g , £ = £ H , p = (e — l)v/c 

z j LO 



Under a Lorentz transformation, with accuracy to the term linear in v/c, we 
have u/ = uj — q x v and q' = q— xuv/c 2 . Note also that p is invariant under 
the Lorentz transformation, i.e. p = p' . The last term in (}2"4*j) gives rise to 
a coupling between s— and p— polarized waves. However, it can be shown 
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18j that this coupling gives a corrections ~ {v/c) 2 to the frictional drag force 
between the layers, so this term can be omitted. The solution of the Maxwell 
equations in the K' reference frame can be written as 



E'= 



Ve ipz + w'e- ipz , 



u 2 e ipz , 



< z < d 
z > d 



(25) 



^From the boundary conditions for the s— and p— polarized waves, which 
follow from we get the equations 



11! 



2ipd 7 / 



■ #2 P (qW) e 

R 2s ^,J)e 2ipd v' n , 



ipd 



47Tpj/ 2 q/e 

£a/(e 2p + 1) 

4Wj /2n ,e^ 
pc 2 (e 2s + 1) 



(26) 
(27) 



The relations between the fields in the K and K' reference frames are deter- 
mined by the Lorentz transformation. As it was shown in Ref. |18j . such a 
Lorentz transformation gives terms of the order v/c which couple the s— and 
p— polarized waves but this result in a contribution to the frictional drag of 
the order (v/c) 2 . Thus we can take this transformation in zero order in v/c 
so that v',(uj') = Vq(u), v' n ,(u') = (uj'/u) v n (u) and similar equations for w. 
After the transformation the solution of the system of the equations (|2()| l2~T| 
l2*7j) take the form 



(28) 

(29) 

(30) 

(31) 
(32) 



W r 



Airp 


j /2q , (q',u/) &P d R lp (q,u) 


J fin (q,w) 


A p 




(e 2p (q',o/) + l)u/ 


(e lp (q,w) + 


Airp 


J/lq (q,uj)e 2i P d R 2p (gW) 


j /2q , (<?>') e ipd 


A p 




(e lp (q,w) + 1)lo 


(e 2p (q',a/) + l)a/ 


4ttuj 


~]f2n> {q'^')^ d R ls {q,u 


) J/in (g, w) 


A s pc 2 


[ (e 2s (q',u') + l) 


(e is (q,u) + 1)_ 


4:7100 




j/io(q,w) e 2ipd R 2s (?W 


) ^(q'^Oe^ 


A s pc 


2 


(e ls (q',a;) + l) 


(e 2s (q',u/) + l) 






qvq 








V z = W z = 






P 


P 



where we have introduce the notation 



A, 



1 - e 21 ^ (qV) R lp (q,u) 
l-c 2ipd R 2s (q',u/) i*i a (q,w) 
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3 Calculation of the frictional drag force be- 
tween 2D systems 

The frictional drag stress a which acts on the conduction electrons in layer 1 
can be obtained from the xz— component of the Maxwell stress tensor a^, 
evaluated at z = ±0 

1 r +oc r i 

Here the (...) denote statistical averaging over the fluctuating current densi- 
ties. The averaging is carrying out with the aid of (jl5j) for s— polarized waves 
and the similar equation for p— polarized waves. Note that the components 
of the fluctuating current density j fi and j j 2 refer to different layers, and are 
statistically independent, so that the average of their product is zero. Ex- 
panding the electric field and magnetic induction in Fourier series we obtain 

° = ^/^{[e<^(q,w)^(q,w)> + (S,(q >W )B:(q 3£ «;)> 

07T J (271") 

+c-4 =+0 -[...U_o} (34) 

For a given value of q it is convenient to express the component E x and B x 
in terms of the components along the vectors q and n 

E x = (q x /q)E <l -(q y /q)E a (35) 
B x = (q x /q)B <l -(q y /q)B n (36) 



After substitution of expressions (j35H36|) into and taking into account 
that the term which is proportional to q y is equal to zero [TH], we obtain 

+cc] 2=+0 - [...U_ } (37) 

where 

E z (z = +0) = (v z + w z ) = (q/p)(w n - u q ) = (qp*/ \ p | 2 )(w q - u q ) (38) 

E z (z = -0) = u lz = {q/p)uq = {q/p){w n + u q ) (39) 
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E^z = +0) = E^z = -0) = V(i + 



w r 



B z (z = +0) = (qc/u)(v n + w n ) = B z (z = -0) = {qc/u)u Xn 
Bq(z = +0) = (pc/u)(w n - v n ) 
B^{z = -0) = (pc/u)u ln 
After substituting these expressions into formula (J37|) we obtain 



(40) 
(41) 
(42) 
(43) 



1 



a 



167T 3 JO 



+oo 



du; / d 2 qq x 



P 



(P + P*)«l^q| 2 >-<K| 2 > 



(I ^q + ^q | 2 )) + "V5)) 

^(p + P'XflWn | 2 )"(hn | 2 )-(k n + Wn | 2 )) 



+ 1- 



-(p-p*)((Un<-UnO)]) 



(44) 



where we integrate only over positive values of tu, which gives an extra factor 
of two. 

Substituting (|2~%j) and (|3^j) into (j4^j) and taking into account that p = p* 
for q < ijjjc and p = —p* for q > u/c, we obtain 



ft 

?7T 3 JO 



. d 2 gg a 



x 



<j<(oj/c)v / £ 

ri p (^)T 2p (u; - q x v)(n(u - q x v) - n{u)) 
| 1 - e 2i P d R lp (uj)R 2p (uj - q x v) | 2 



T 1p (uj)(\ 1-R lp (u) 


2) + (| l-e%( W ) 


2 )(nH + l/2) 




1 - e 2i P d J R lp (w) J R 2p (cj - q x v) 


I 2 



2?r 3 Jo 



dcu 



q>(w/c)y/e 



-2|pM 



X- 



Im_R lp (co')Im_R2p( ( ^ — 9a: 

1 -e^bl^H^p^-g^ 



-(n(w - q x v) -n(u)) 



+ {P 



(45) 



where 



T ip (uj) 


= 1- 


1 R I 2 

-T^ip \ 


1 1 


p 12 
-'Hp I 




= 1- 


1 R I 2 


1 


p \2 



167rReo"z(a;)p 
a;£|eiz + 1| 2 

167rRecJi(a;)a; 
pc 2 \e it + 1| 2 
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The first integral in (|66|) is the contribution to the frictional drag force from 
propagating electromagnetic waves. This integral contains terms which for- 
mally diverge upon integration over u. These terms are proportion to uj^ 1 
at large frequencies and appear as a result of the expansion of the reflection 
factor R 2 p( s )(uj — q x v) in power series in q x v, and upon performing the q in- 
tegration. A similar divergence also occur in the derivation of the static van 
der Waals interaction [T7] . and result from zero point vacuum fluctuations 
of electromagnetic field. The solution of this problem consists of subtraction 
from the integrand the terms which do not depend from separation between 
layers d in the limit d — > oo. In our case this procedure consists of subtraction 
from the first integrand in (|66|) the same expression taken at T = K and 
with denominator equal to unity. The second term in (|66|) is derived from 
the evanescent field. 



4 Some limiting cases 

Consider distances d « dw ~ ch/ksT (at T = 3 K we have d\y ~ 10 6 A). 
In this case we can neglect by the first integral in ()66|) . put p ~ ig and extend 
the integral over q to the whole q— plane. Using these approximations, the 
second integral in ()66|) can be written as 



a 



h 

X 
X 



+oo 



dq y / dq x q x e 2qd 

» JO 

o 

du[n(u) — n(uj + q x v)] 

lmRi p (u)lmR2 P (uj + q x v) 
1 _ e-^ d R lp {u)R2 P {u + q x v) | 2 



P) 



+ s 



duo[n{uo) + 1/2] 
P)]} 



lmRi p (uj - q x v)lmR 2 p{uj) 
1 - e-^ d R lp {u - q x v)R 2p {u) 



(46) 



The second term in this expression is proportional to v 2 as v — * and can be 
neglected in the limit of small v. In the first term we can use approximation 



n(u) - n(u + q x v) -q x v— 



dn e Ku/k B T hqxV 
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Thus 



a 



Tiv 
2^ 



dgg 3 e~ 




ImRip (uj) Im_R 2p iyj) 
1 _ e-WR lp (u)R2 P (u) 



+ [P 



Note that a is linear in the velocity v. Let us describe the 2D layers in RPA 
approximation. Thus for q < k F (corresponding to separations d > kp 1 , where 
kp is the Fermi wave vector of the degenerate electron gas system; for 2D 
electron layer with electron density n s « 1.5 • 10 n cm~ 2 , kp = (27m s Y^ 2 ~ 
10 6 cm _1 ) the transverse and longitudinal parts of the conductivity for 2D 
electron layer can be written in the form 



itoe 2 n s { uju 1 

° l = 7 . : r-r—r • 1 48 

q 2 e F { (u) + 17) V u 2 - 1 - ryu J 

2ie 2 n s U(v / H 2 - 1 - H) .... 

0t = 71 — r^i ( 49 ) 

m*[uj + 17) 

where U = (uj + vy)/qvp, 7 = 1/r, Vp = hk F /m* is the Fermi velocity, r is 
a relaxation time, e F = Ti 2 k 2 F /2m* is the Fermi energy. In the experiment 
0| m* = 0.067m e ,v F = 1.6 x 10 7 cm/s, ep ~ 60 K and the mobility jj, ~ 
2xl0 6 cm 2 /Vs, so that r ~ 7.6 x lCT n s. Let us divide the integration over 
< q < 00 into the two parts < q < uj/vp and cu/vp < q < 00. In the 
first part of integration u > 1, and taking the limit m >> 1 we obtain in this 
limit the Drude formula for conductivity 

ie 2 n s 

= <?t = — -, r (50) 

m*(uj + 17) 

In the second part of integration u < 1 and taking the limit u « 1 we 
obtain 

^ = -5— (51) 



e 2 n s vp 



e F q 

where we put 7 equal to zero because it gives only a small contribution in 
this limit. 
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Let us consider the case of small separation d when a = (2kBTd/hvp) < 1. 
Introducing the dimensionless variables q = x/2d and u = (kBT/h)y we 
obtain in this limit for ay < x < oo 

M* + i°l*) w i- - (52) 

ar + A p (x + iay J A p (x + lay) 

R, = (53) 
iA s y - x l 



and for < x < ay 



X ' P X 

Rp = X' p x-2y 2 -2iy5 {M) 
3 2xy + \' s y + i2x6 1 j 



where 



8ne 2 n s d ( e 2 n s \ , 27m s e 2 ( H \ 2 , 8nn s e 2 d 

em Vp \m*c z J em*d \Kb1 J m c A 

We note that the expression (J53)l has pole at 

2 2irn. q e 2 



u 2 = & —q (56) 

em* 

what corresponds to the plasmon excitations [2E|- After substituting (|52T - l5l)j) 
in (}4Tj) we obtain for the frictional drag rate 

1 (kT) 2 fk B T\ 5 (k B T\ 2 . . 

where and are the contributions from s— and p— polarized waves, 
respectively, Qtf — 2e 2 m* /h 2 e is the single-layer Thomas- Fermi screening 
wavevector, €tf = h 2 q^ F /2m* .The first term in (|57j) agrees with the result 
of Gramila at a/ [3] and Persson and Zhang [2Z1- From comparison ()57|) and 
(jHHJl it follows that for 
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the contribution from p— polarized waves exceeds the contribution from s— 
polarized waves for all distances d < hvp/ksT. However for n s > n c the 
contribution from s— polarized waves will dominate for d > lQ(e/n s ) l l 2 . 
For example, for T = 3 K and for the conditions of experiment of Ref. [3 E] 
n c ~ 10 12 cm~ 2 and we find that in this case the retardation effects are 
small. However retardation effects are important for high electron densities. 
For example, assuming that e — 1 and n s Pd 10 15 cm~ 2 , which correspond to 
about 1 monolayer of silver, we find that the contribution to frictional drag 
from the retardation effects will dominate for d > 15 A. This is illustrated in 
Fig. 2a which shows the shear stress when the relative velocity v = lm/s. We 
have performed calculations for two different temperatures, T = 273K and 
77K, and the s and p-wave contributions are shown separately. In Fig. 2b 
we show the same quantity for two quantum wells at T = 3K and with n s = 
1.5 x 10 n cm~ 2 , m* = 0.067m e , vp = 1.6 x 10 7 cm/s and r = 7.6 x lCT n s, and 
with e = 1. In this case the p-wave contribution dominates for d < 1000 A. 

Let us estimate the voltage JJ\ induced in a thin silver film (layer 1) 
(open circuit) when a current flow in another parallel silver film (layer 2). A 
voltage difference of order 1 pV can be measured with standard equipment so 
that if Ui is of order pV or larger, it is possible to probe retardation effects 
with this experimental setup. If L denote the length of the metallic films 
(assumed identical) in the direction of the driving current, then U\ — LE\ 
and XJi = LE 2 = LJijoi where o 2 = n 2 e 2 r 2 /m* is the conductivity (r 2 is 
a Drude relaxation time and m* the electron effective mass). Thus, using 
the equation (see introduction) 7 = nin2e 2 Ei/J 2 with E\/J 2 = Uxjip^'i) = 
(Ui I U2)m* I (n 2 e 2 T 2 ) gives U\ = (jT/rrfnx)^- In a typical case r = 4 x 
10 -14 s and n\ ~ 10 15 cm -2 , and from Fig. 2a, 7 ~ 10~ 6 Ns/m 2 , giving 
Ui ~ 10 -8 U 2 . Thus if the applied voltage U 2 ~ 1 V, the induced voltage 
would be of order 10 pV, which should be possible to measure. 



5 Frictional drag between 3D systems 

For high electron densities, when the thickness of the layers h >> -nr 1 / 3 , 
where n is a volume electron density, the electrons behave as in 3D systems. 
It was shown in Ref . jTHj that for 3D systems the frictional drag stress is also 
given by formula (|47|). where the electromagnetic reflection coefficients 

R ip = ■ , R is = — ■ — , (60) 

Etf + esi p+ Si 
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where £j is the complex dielectric constant for layer i, 

[up 

s * = y^-<7 2 ( 61 ) 

Consider two identical 3D layers described by the dielectric function 



uj 2 



= 1 - , *. - y (62) 

where r is the Drude relaxation time and u p the plasma frequency. For 
k B T > h'j, for small frequencies and d < c/to p (h'f /ksT) 1 ^ 2 



ImRp « — , ReRp 

U 2 T 

^ ^ t t~~\ fUp\ 2 u) „ 9 / co>„ \ 2 u; 

ReR s « 0, Imi? s w4 M — r, for g 2 > — -, 

V c / 7g 2 V c / 7 

then, taking into account that for 3D systems n s = nh, (jlTj) gives 



2 



.-1 



- _ 13 32 ( £ ^ T ) 

Dp he F {k F d) 2 {k TF d) 2 {k F h){u p r) 2 1 J 

_ e\k B T) 2 M 2 

where A;f._p = 6irne 2 /e F is the 3D Thomas- Fermi screening wavevector, and 
k F = (37r 2 ra) 1//3 is the 3D Fermi wavevector. From comparison (JB^j) and 
(JMj) it follows that, for d > £ 1 / 2 (c/ 'up)(u p T)~ 1 , the s— wave contribution 
exceeds the p— wave contribution Thus, in the case of the dissipative van der 
Waals interaction between 3D bodies, retardation effects become important 
for much shorter distances in comparison with conservative one, when the 
retardation effects become important for d > c/u p |17j . 

Fig. 3 shows the calculated shear stress for two semi-infinite silver bodies 
moving with the relative velocity v = lm/s parallel to the flat surfaces. 
Results are shown for the s- and p-wave contribution, where in the latter case 
we have taken into account non-local effects [the dashed lines shows the result 
when the local (long-wavelength) dielectric function is used]. Results are 
shown for two different temperatures, T = 70K and 300K, and the observed 
temperature dependence reflect both that of the temperature prefactor T 2 in 
the expression for the shear stress, as well as the temperature dependence of 
the 
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6 Relation between friction and heat transfer 

The frictional shear stress studied above is closely related to the heat transfer 
from one solid to another when the solids have different temperatures. For 
large separation the heat transfer is given by the Stefan's law 

j 2 = ^L( Tl 4 -T 2 4 ) (65) 

where T\ and T 2 are the temperatures of solid 1 and 2, respectively. This 
formula correspond to emission of real photons. However, for short separation 
d it is possible for the evanescent near field to transfer energy from one solid 
to the other. This correspond to photon-tunneling. In general, the heat flux 
(energy flow per unit area and unit time) is given by a formula very similar 
to that for the frictional stress fl§\ I3U] : 



poo p 

J z = — / duju / d 2 q 

87T J JO Jq<uj/c 



(1- 




2 )(1- 


R 2p (uj) 


2 )(ni(u;) - n 2 (uj) 




1 


-e^ d R lp (u)R 2p (u) 


2 



— f duu [ d 2 qe~ 

ZTT JO Jq>uj/c 

lmRi p (u) Im_R 2 r, (uj) 



■2\p\d 



2p[ 



I 1 _ e- 2 \p\ d R lp (u)R 2p (u) | 2 
+ [p -> s] 



lni(u) - n 2 (uj)) 



(66) 



where 



hwjk B Tx 



1 



(67) 



is the Bose-Einstein factor of solid 1 and similar for n 2 . Fig. 4a shows 
the heat transfer between two semi-infinite silver bodies separated by the 
distance d and at the temperatures T\ = 273K and T 2 = OK. The s and 
p-wave contribution are shown separately, and the p-wave contribution has 
been calculated using non-local optics (the dashed line shows the result using 
local optics). It is remarkable how important the s-contribution is even for 
short distances. The detailed distance dependence of J z has been studied by 
Ploder and Van Hove within the local optics approximation, and will not be 
repeated here. The nonlocal optics contribution to (J z ) p , which is important 
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only for d < I (where / is the electron mean free path in the bulk), is easy to 
calculate for free electron like metals. The non-local contribution to IhlR p is 
given by [27] 

Using this expression for IhlR p in (jUol ) gives the (surface) contribution: 
where 

POO / 1 1 \ 

= 0.1827 / dyy 3 (- j^y- -> 1.186 

as T2/T1 — > 0. Note from Fig. 4a that the local optics contribution to (J z ) p 
depends nearly linearly on 1/d in the studied distance interval, and that 
this contribution is much smaller than the s-wave contribution. Both these 
observations differ from Ref. [31], where it is stated that the s contribution 
can be neglected for small distances and that the p-wave contribution (within 
local optics) is proportional to 1/d 2 for small distances. However, for the 
very high-resistivity materials, the p-wave contribution becomes much more 
important, and a crossover to a 1/d 2 - dependence of (J Z ) P is observed at very 
short separations d. This is illustrated in Fig. 4b and 4c, which have been 
calculated with the same parameters as in Fig. 4a, except that the electron 
mean free path has been reduced from I = 560 A (the electron mean free 
path for silver at room temperature) to 20 A (roughly the electron mean 
free path in lead at room temperature) (Fig. 4b) and 3.4 A (of order the 
lattice constant, representing the minimal possible mean free path) (Fig. 4c). 
Note that when I decreases, the p contribution to the heat transfer increases 
while the s contribution decreases. Since the mean free path cannot be 
much smaller than the lattice constant, the result in Fig. 4c represent the 
largest possible p-wave contribution for normal metals. However, the p-wave 
contribution may be even larger for other materials, e.g., semimetals, with 
lower carrier concentration than in normal metals. This fact has already 



18 



been pointed out by Pendry: the p-wave contribution for short distances is 
expected to be maximal when the function 



ImRp ~ 



Im- 



e+1 



Im 



2— — + 



UJpT / 



is maximal with respect to variations in 1/r. This gives: 

2k B T 



LdpT = 



flUOr 



where we have used that typical frequencies uo ~ ksT/h. Since the DC 
resistivity p = A-n/iuJ^r) we get (at room temperature) p ~ 2nh/kBT ps 
0.14 Ocm. 



7 Summary and conclusion 

We have used a general theory of a fluctuating electromagnetic field to cal- 
culate the frictional drag force between 2D and 3D electron systems. The 
separation d between the parallel electron layers is assumed to be so large 
that the only interaction between the layers is via the electromagnetic field 
associated with thermal and quantum fluctuations in the layers; the resulting 
friction force can be considered as the dissipative part of the van der Waals 
interaction. A general formula has been obtained, in which the frictional 
drag force is expressed through the electromagnetic reflection coefficients for 
s and p— waves. We have found that the non-retarded Coulomb interaction, 
connected with evanescent p— polarized waves, is the predominant process 
for small layer separations and small electron densities. For high electron 
densities retardation effects (connected with evanescent s— polarized waves) 
become very important, and we have suggested a new experiment, involving 
thin metallic films, where the theory can be tested. We have shown that 
retardation effects are even more important for interaction between 3D elec- 
tron systems. For very large separations the interaction is dominated by 
the traveling electromagnetic waves, which results from black-body radia- 
tion. However, the latter interaction appears negligible in comparison with 
phonon mediated process. Finally, we have pointed out the close relation 
between heat transfer and friction. 
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FIGURE CAPTIONS 

Fig. 1 Left: a metallic block sliding relative to the metallic substrate 
with the velocity v. An electronic frictional shear stress a will act on the 
block (and on the substrate). Right: The shear stress a can be measured 
if instead of sliding the upper block, a voltage U 2 is applied to the block 
resulting in a drift motion of the conduction electrons (velocity v). The 
resulting frictional stress on the substrate electrons will generate a voltage 
difference U\ (proportional to a) as indicated in the figure, which can be 
measured experimentally. 

Fig. 2 The shear stress as a function of the distance d between the surfaces 
(the log-function is with basis 10). (a) For ^monolayer- films of silver for two 
different temperatures. The s and p-wave contributions are shown separate. 
In the calculation r = 4 x 10~ 14 s and 20 x 10" 14 s for T = 273K and 77 
K, respectively. We have assumed n s = 1.05 x 10 19 m~ 2 , m* = m e , and 
vp = 1.4 x 10 6 m/s. (b) For quantum wells at T = 3 K. In the calculation 
t = 7.6 x 10~ n s, n s = 1.5 x 10 15 m" 2 , m* = 0.067m e , and v F = 1.6 x 10 5 m/s. 

Fig. 3 The shear stress as a function of the distance d between the surfaces 
of two semi-infinite silver bodies. The s and p wave contributions are shown 
separately and for two different temperatures, T = 70 K and 300 K. The 
p-wave contribution has been calculated both using a local dielectric function 
(dashed lines) and using a theory which takes into account nonlocality within 
the jellium model. 

Fig. 4 (a) The heat transfer flux between two semi-infinite silver bodies, 
one at temperature 7\ = 273 K and another at T 2 = K. (b) The same as 
(a) except that we have reduced the electron relaxation time r for solid 1 
from a value corresponding to a mean free path vpT = I = 560 A to 20 A. 
(c) The same as (a) except that we have reduced / to 3.4 A. 
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